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Event shape data from e"'"e~ annihilation into hadrons coUected by the JADE experiment 
between -^s = 14 and 44 GeV are used to determine the strong coupling as. QCD predic- 
tions complete to next-to-next-to-leading order (NLLO), alternatively combined with next-to- 
leading-log-approximation (NLLA) are used. The stability of the NNLO and NNLO+NLLA 
results with respect to variations of the renormalisation scale is improved compared to pre- 
vious results obtained with next-to-leading-order (NLO) or NLO-f-NLLA predictions. The 
energy dependence of as agrees with the QCD prediction of asymptotic freedom and excludes 
absence of running with 99% confidence level. 



1 Introduction 

The determination of the stron g couphng as in Quantumchromo Dynamics (QCD) is an im- 
portant test of the theory Hl^^^. The analysis of hadron production in e'''e~ annihilation with 
observables based on jet production and event shape definitions plays a major role, because in 
the e^e~ environment there is no interference between intial and final state and one generally 
has a clean experimental environment with small background rates and neglegible pile up. Re- 
cently the NNLO corrections for QCD predictions of event shape observable distributions'^'^ and 
moments^ were completed and used for a measurement of ag using data of the LEP experiment 
alephCI. It is therefore of high interest to measure as with JADE data and the same improved 
QCD predictions'^. We also include matching of existing NLLA QCD calculations with the new 
NNLO calculations^]. 

The JADE experiment operated at the PETRA e+e" colhder at DESY from 1979 to 1986 
at centre-of-mass (cms) energies ^/s = 12 to 44 GeV. Significant data samples were recorded 
at ^ = 14,22,34.6,35,38.3 and 43.8 GeV. The JADE detector is described in^O. The JADE 
experiment combines tracking of charged particles in a solenoidal magnetic field and electro- 
magnetic calorimeters over a large solid angle. This allows to select hadronic final states with 
high efficiency while demanding that the events are fully contained in the detector. The JADE 
reconstruction, simulation and analysis software was ported to current computers^] and is used 
in our work. 

At the JADE energies background processes relevant for hadron production are e^e~ — 
r+r^ with subsequent hadronic r decays and two-photon interactions with hadronic final states. 
Initial state radiation (ISR) of photons from the beam particles also occurs with a significant 
rate leading the production of hadronic final states at reduced cms energy. After application 
of suitable event selection cuts based on 4-momentum balance and particle multiplicity the 
background processes can be neglegted. Residual effects of ISR are corrected for. 

We use the event shape observables thrust 1 — T, heavy jet mass Mh, total and wide 



broadening and i?W) C-parameter C and 1/23, the value of the distance definition yij = 
2 min(£'j^, Ej){l — cos 6ij)/s of the Durham jet algorithm where the event changes from a three- 
to a two-jet configuration^l2lSl_ 

2 Data Analysis 

The distributions of event shape observables are computed from all selected events, and from 
samples of simulated events at several levels. The final state after termination of the parton 
shower is called parton-level, the final state after hadronisation and decays of particles with life- 
times shorter than 300 ps is called hadron-level, and after passing the events through the JADE 
detector simulation, recontruction and analysis procedures the result is called detector-level. We 
use the programs PYTHIA, HERWIG and ARIADNE tuned by the OPAL collaboration for this 
purpose, see^ for details. 

The contribution of e'''e~ — > bb events is subtracted from the distributions using simu- 
lated events at detector-level. The results are multiplied by the ratio of simulated hadron- and 
detector-level distributions to correct for experimental acceptance and resolution. In this way 
we obtain distributions of event shape observables corrected to the hadron-level. A compari- 
son of observed and corrected event shape distributions with predictions by the Monte Carlo 
simulations shows satisfactory agreement. 

3 Determination of ag 

The hadron-level distributions are compared with the QCD predictions corrected for hadroni- 
sation effects by a x^-method with as as a free parameter and the full statistical covariance 
matrix. The hadronisation corrections for the QCD predictions are found for cumulative the- 
ory predictions by dividing cumulative distributions at hadron- and parton-level. In order to 
justify this procedure we compare simulated distributions at parton-level with the NNLO QCD 
calculations for as(-^zo) = 0.118 and find reasonable agreement as shown in figured] (left). The 
program PYTHIA is used in the standard analysis and HERWIG and ARIADNE are used as 
alternatives. 
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Figure 1: (left) Ratio of QCD and MC prediction at parton-level for observables and cms energies as indicated 
on the figures. The different line types show the ratios for different MC programs, (right) Data at hadron-level 
for 1 — T with superimposed NNLO-f NLLA fits at the six JADE energy points. The arrows correspond to the fit 

ranges. 



Only restricted regions of the distributions are used for the fits, i.e. those where the ex- 
perimental and hadronisation corrections are stable and where the QCD predictions have LO 
contributions and don't have large logarithmically enhanced terms The resulting fit ranges 
are the same for all energy points and are shown on figure [T] (right) for 1 — T. 

The experimental uncertainties are found by repeating the analysis with varied selection 
cuts and procedures. The uncertainties from the hadronisation corrections are determined by 
comparing fit results obtained with PYTHIA, HERWIG or ARIADNE. The theoretical uncer- 
tainties are found by changing the renormalisation scale parameter = [ij \fs from its standard 
value = 1 to = 0.5 or 2.0. 

Figure [T] (right) shows the NNLO+NLLA QCD predictions using the fit results for ag su- 
perimposed over the hadron-level data for 1 — T. The agreement between data and theory is 
excellent for both NNLO and NNLO+NLLA fits. For example, the x^/d.o.f. values for 1— T with 
NNLO fits range from 0.7 (14 GeV) to 2.5 (34.6 GeV). The root-mean-square (RMS) values of 
as results from NNLO+NLLA fits to the six observables are between 0.007 (22 GeV) and 0.003 
(44 GeV), consistent with the theory uncertainties of the combined values (see below). This is 
an important cross check made possible by having six different observables ii i the analysis. 

We combine the results from the six observables at each energy point as in^lSISl ^nd show the 
results of the NNLO analysis in figure [2] together with the results o The result of combining 
the combined values of as after evolving to the reference scale M^o is shown as lines in the same 
figure. The final result using this combination is for NNLO 

as(Mzo) = 0.1210 ± 0.0007(stat.) ± 0.0021(exp.) ± 0.0044(had.) ± 0.0036(theo.) 

and for NNLO+NLLA 

as(Mzo) = 0.1172 ± 0.0006(stat.) ± 0.0020(exp.) ± 0.0035(had.) ± 0.0030(theo.) . 




Figure 2: (left) Combined values of as from NNLO fits at the JADE energy points together with the results o 
The lines show the evolution of as based on the combined result, (right) Combined values of as(M2o) for each 
analysis. The shaded bands and dashed lines show the combined values of as(M2o) with total uncertainties. 

Our final result from NNLO+NLLA predictions has the smaller uncertainties of 4%. It is 
consistent with previ ous d eterminations of a<s,{Mi<3) using NNLO theoryl^and with with recent 
averages of as(M2o)^I^2EI_ Repeating the combination of the combined NNLO+NLLA results 
at each energy point without evolving to a common reference scale results in a probability of 



10~^ assuming partially correlated hadronisation uncertainties and ignoring theory uncertainties. 
We interpret this as strong evidence for the dependence of the strong coupling on cms energy 
as predicted by QCD from JADE data alone. 

Comparison with NLO and NLO+NLLA fits using the same data, hadronisation corrections 
and fit ranges gives consistent results. The theory uncertainty of the NLO+NLLA analysis is 
larger by about 60% compared with our NNLO or NNLO+NLLA results. The theory uncertainty 
of the NLO analysis with fixed renormalisation scale = 1 is larger by a factor of 2.6 compared 
with the NNLO analysis. Figure [2] (right) shows a summary of the different analyses. The 
smaller theory uncertainties of the NNLO or NNLO+NLLA analyses compared with NLO or 
NLO+NLLA analyses follow from the weaker dependence of the NNLO calculations on the 
renormalisation scale parameter x^^. 

4 Summary 

We have shown a determination of the strong coupling constant 03(^20) with NNLO and 
NNLO+NLLA QCD calculations using JADE data. The final result as(Afzo) = 0.1172 + 0.0051 
has a competetive uncertainty of 4% and is consistent with other measurements. Our data 
confirm the running of as as predicted by QCD. 
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